Fruit flies modulate passive wing pitching to generate in-flight turns 
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Flying insects execute aerial maneuvers through subtle manipulations of their wing motions. Here, 
we measure the free flight kinematics of fruit flies and determine how they modulate their wing 
pitching to induce sharp turns. By analyzing the torques these insects exert to pitch their wings, 
we infer that the wing hinge acts as a torsional spring that passively resists the wing's tendency 
to flip in response to aerodynamic and inert ial forces. To turn, the insects asymmetrically change 
the spring rest angles to generate rowing motions of their wings. Thus, insects can generate these 
maneuvers using only a slight active actuation that biases their wing motion. 



To generate the vertical force necessary to sustain 
flight, small insects must beat their wings hundreds of 
times per second. Under this constraint, how do they 
manipulate these fast wing motions to induce flight ma- 
neuvers? Although recent studies have made progress ad- 
dressing how wing motions generate aerodynamic forces 
[H 0- 3, Si , understanding how the wing motions them- 
selves arise and what control variables govern them re- 
mains a challenge. To address these questions, we ana- 
lyze the torques freely-flying fruit flies (D. melanog aster) 
exert to move their wings. Specifically, we elicit sharp 
free-flight turns from these flies and measure their wing 
and body kinematics. By using a model of the aerody- 
namic forces on flapping wings, we extract the torques 
the insects exert to generate the wing motions. From 
these torques, we construct a mechanical model of the 
wing rotation joints that demonstrates how the interplay 
of aerodynamic, inertial and biomechanical forces gener- 
ate the wing kinematics. Finally, we connect this model 
to the turning dynamics of flies and describe the wing 
actuation mechanism that unifies these maneuvers. 

To quantify turning kinematics in fruit flies, we first 
use three orthogonal cameras to capture their free flight 
maneuvers at 8000 frames per second or about 35 frames 
for each wing beat. We then reconstruct the three- 
dimensional wing and body motion of the flies from these 
videos using the motion tracking techniques described in 
The body kinematics, described by the centroid co- 
ordinates and three Euler angles - yaw, body pitch, 
#6, and roll, are shown in Fig. [T]A. and visualized in 
Fig. [TJ>. During the level flight, the fly performs a 120° 
turn in 80 ms, or 18 wing beats. 

To induce such a turn, the insect generates differences 
between the motion of its left and right wings. We quan- 
tify these changes by plotting in Fig. [T]G the time course 
of three Euler angles - stroke, </>, deviation, 0, and wing 
pitch, ip - that describe the orientation of the wings rela- 
tive to the hinges they rotate about 0. A three dimen- 
sional representation of a typical wing stroke is shown 
in Fig. [Tp. Remarkably subtle asymmetries between the 
pitch angles of the wings drive the turn. For a discus- 



sion of why the other observed asymmetries produce a 
negligible effect on the yaw dynamics see [7]. To quan- 
tify the wing pitch asymmetry that induces the maneuver 
we phase- average [8[ six consecutive strokes from t = 10 
ms to 30 ms, during which these asymmetries are most 
prominent. We find that the symmetric front- and back- 
strokes of the left wing result in canceling drag forces 
and no net torque on the fly. In contrast, the right wing 
beats with mid-stroke angles of attack a = 49° during 
the front- and a = 40° during the back- stroke. The 9° 
difference in a results in a net drag force towards the 
back of the fly and a torque that turns the fly clockwise. 
Thus, by varying wing pitch and consequently angle of 
attack, insects row through the air to perform a turn. 

To elucidate how the fly controls its wing pitch, we 
invert the equations of motion of the wings to determine 
the torque the fly exerts at the wing hinge, 



• uj — r x m w a 



(i) 



From kinematic data, we determine the rotational accel- 
eration Co and acceleration a of the wing centroid. We 
also measure the wing mass m w and use morphological 
measurements of fruit fly wings to estimate the wing mo- 
ment of inertia I w and center of mass to hinge vector r 
[9| . The aerodynamic torque about the wing hinge, r a , is 
calculated using a quasi-steady model detailed in [7|, LL0[ . 
From Eq. (pp), we determine the torque component flies 
exert to pitch their wings, r p , during the turn and during 
steady flight. 

Previous studies have suggested that the pitching 
torques result from torsional deformation of the wing 
hinge and simply resist the tendency of aerodynamic and 
inertial forces to flip the wing 0, [n|, E3, Gj, Q Est 
To integrate these ideas into a dynamical model for the 
pitching torque, we plot r p versus i/j in Fig. [2]A. for the 
final 9 consecutive strokes associated with steady flight 
in Fig. [T]0. The torque data traces out an elliptical curve 
whose major axis has a negative slope. This negative cor- 
relation indicates that when the wing angle deviates from 
approximately 90°, the hinge produces a restoring torque. 
The area enclosed by the ellipse indicates the energy dis- 
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FIG. 1: (A) Measured fruit fly body position and orientation versus time. (B) Phases of motion of a fruit fly performing a 
120° turn. The three panels show 6 of 821 frames recorded by high speed videography, while the corresponding reconstructed 
kinematics are visualized using the model flies. Labeled on the fifth phase are the Euler angles that quantify the fly's orientation. 
(C) Measured fruit fly wing orientation relative to body versus time. (D) The trajectory and orientation of wing chords during 
flapping stroke. The red and blue lines on the globe depict the chord's orientation at equal time snapshots, with the ball 
depicting the wings leading edge. Labeled on the globe are the Euler angles used to quantify wing orientation. (E) The 
unwrapped ball-and-stick diagram highlights the 9° asymmetry in the mid-stroke angle of attack that causes the maneuver. 



sipated by the hinge as it pitches the wing. These data 
suggest the hinge acts like a damped torsional spring, 



-ft (^ - ijj ) - Cip, 



(2) 



where the parameters ft, C, and correspond, respec- 
tively, to the torsion constant, damping constant and the 
rest angle of the torsional spring. We fit this model to 
all 9 wing strokes, finding values ft = 91 ± 9 pNm/°, 
C = 39 ± 12 fN m s/°, and i/j = 90 ± 1° that account for 
95% in the variance of the pitching torque (Fig. [2]/Y). 

To determine how the wing pitch is actuated differently 
during the turn, we repeat the above analysis for the 6 
strokes that induce the maneuver. We find that Eq. ([2]) 
also accounts for the exerted torques during these wing 
beats. We compare the phase- averaged torque for the 6 
asymmetric strokes with the phase- averaged torque for 
the 9 symmetric steady strokes in Fig. [2j3. We find that 
the two loops are shifted horizontally with respect to each 
other, indicating a change in ^o, the wing's rest angle. 
In fact, by plotting the values for ft, C and as a func- 
tion of time, we show that only tpo varies significantly 
throughout the maneuver (Fig.[3]A.-C). The steadiness of 
ft and C suggests that these parameters correspond to 



material properties of the wing hinge. This interpreta- 
tion is supported by the fact that the value of the spring 
constant, ft, agrees with scaled estimates for the torsion 
constants taken from measurements on the wings of house 
flies [12|. Comparison of the parameter values with the 
fly's yaw dynamics in Fig. [3p indicates that, to induce 
the clockwise turning maneuver (orange shading), the in- 
sect increases the of the right wing relative to the left 
by about 15° for 6 strokes. We find that for the two sub- 
sequent wing beats (blue shading) the insect decreases 
the ipo of the right wing relative to the left by about 10°. 
This reversal in the sign of A^o = — Wf^ indicates 
that the insect is generating a counter-clockwise torque 
that slows its yaw velocity. In the final 9 wing strokes 
(gray shading), the values for the left and right wings 
are nearly equal so that no active torque is generated. 

To validate that the yaw dynamics of the fly derives 
from changes to the relative rest angle of the wings, 
A^o, we simulate the coupled wing-body dynamics in 
a model fly. The driving of each wing is simulated by 
prescribing its stroke and deviation angles. We then de- 
termine the fly's yaw angle and wing pitch angles from 
the interaction of the spring model in Eq. ([2|) with the 
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FIG. 2: (A) The pitching torque r p versus the wing pitch ip 
for the wing beats associated with steady flight. Red circles 
represent instantaneous values that correspond to frames of 
the flight sequence. Black triangles are the result of fitting r p 
with a damped torsional spring model. The counter-clockwise 
direction of propagation indicates that the motion is damped. 
(B) The phase- averaged r p versus ip for strokes associated 
with symmetric wing movements (gray trajectory) compared 
to asymmetric wing movements (orange trajectory). 



aerodynamic forces on the wings [7|. To form a mini- 
mal model of the turn, the variables k and C are held 
constant (dashed lines in Figs. [3]A.,B), is set to zero, 
and (j) is driven sinusoidally with flapping amplitude and 
frequency that match the experiments. The rest angle, 
^o, for the left and right wings are prescribed, respec- 
tively, by the blue and red dashed lines in Fig.(3]G, which 
capture the observed trends. Further, we ensure that lo- 
cally the area between the dashed lines matches the area 
enclosed between the two experimental curves. In Fig. 
[3p we compare the observed yaw dynamics of the insect 
to those predicted by our minimal turning model. We 
find that the turn predicted by the model fly (dashed 
line) quantitatively matches experimental measurements 
(solid line) capturing the total accumulated yaw angle 
and time-course of the turn. In Fig. [3p, we compare 
the measured wing pitch with our model predictions (in- 
sets) for the three regions that correspond to A^o < 0, 
A^o > and A^o ~ 0. We find that the minimal model 
recovers the average measured wing pitch asymmetries, 
A^, in these three regions, but does not reproduce the 
detailed sub- wing beat variations in ij;. When we ex- 
tend our simulation to prescribe <p and of the wings 
as measured experimentally, these details are also repro- 
duced [7]. The fact that the turning dynamics of the fly 
are predicted by the minimal model and are insensitive 
to the details of the stroke and deviation angles of the 
wing beat strongly indicates that the turn is controlled by 
Ai/jq. This finding is further underscored by the fact that 
simulation of experimentally measured <p and angles 
without asymmetries in the rest angle, A^o = 0, does 
not lead to asymmetric wing pitch angles and therefore 
does not yield a turn 0. 

Taken together our findings show that as flies flap their 
wings back and forth, the wing hinges resist the tendency 
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FIG. 3: (A)-(C) Torsional spring parameters versus time ex- 
tracted with error bars from the experiment data. (D) Yaw 
angle versus time. The solid-line corresponds to experimental 
measurements and dashed-line to simulation of a model fly 
that flaps with zero deviation angle and a sinusoidal stroke 
angle with amplitude and frequency that match the measured 
values. In the simulation, k and C are held constant as in 
(A) and (B) and ipo is varied as the dashed lines in (C) (E) 
Phase- averaged pitch of the right and left wings, ip, for one 
complete stroke in the three labeled regions. Solid lines are 
experimental measurements, while dashed lines in the insets 
are corresponding strokes from the model fly. The shift Aip 
in the three regions are 6.5° (orange A^o < 0), 2.5° (blue 
A^o > 0), and < 1.0° (gray \Aipo\ < 1°). The corresponding 
Aip for the model fly are 8°, 2.9° and 0°, respectively. 



of aerodynamic and inertial forces to flip the wings. In 
steady flight, this passive mechanism causes the pitching 
of insect wings and results in front- and back- strokes 
with nearly symmetric angles of attack. The rest angles 
of the wings act as control levers that break this symme- 
try (Fig.[3]/Y). This model reconciles how insects combine 
active and passive modulation of their wing kinematics to 
control flight ID, 13, 16|, O ur simulations show 

that the shift in relative rest angles produces a linear shift 
between the wing pitch angles so that Aip « /jl Aipo (Fig. 
HP). We find from experimental analysis of 147 distinct 
wing strokes that \i = 0.51 ± 0.03 which compares well 
with simulations of the sinusoidally flapping model fly 
where /i = 0.6 [7]. The asymmetry in wing pitch angles 
causes asymmetric drag forces on the wings, which are 
related to the fly yaw dynamics by: 



I b 4> b + 20C r <j) b = 2C T 2 A^, (3) 
where I b is the moment of inertia of the fly about the yaw 
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FIG. 4: (A) The minimal model of turning dynamics: the 
wings are driven with a sinusoidal stroke, while yaw and wing 
pitch are determined by aerodynamic, inertial, and spring 
torques. Biasing the rest angle of the spring, produces an 
asymmetric stroke that generates thrust. A difference in rest 
angles, A^o, of the left and right wings induces a turn. (B) 
The average values of Aip versus Aipo for 147 strokes (black 
points). We find Aip = /xA^o, where \i — 0.51 ±0.03 from ex- 
periment, and \i — 0.6 from simulations of the minimal turn- 
ing model. Measurement errors are on the order of the spread 
of the data and are omitted for clarity. (C) The total turn 
angle of the fly versus cumulative pitch asymmetry. Measure- 
ments of ten turns (red dots) are compared with prediction 
of the minimal turning model (dashed line). Measurement 
errors are on the order of the marker size and are omitted for 
clarity. 

axis, C T is a parameter that depends on the wing drag 
coefficient and geometry, uj is the average angular velocity 
of the wings, and fa is the yaw of the body When 
Aip = 0, Eq. (j3]) reduces to the passively damped yaw 
motion characteristic of flapping flight [H, When 
Aifj 7^ 0, the fly generates a driving torque that actively 
yaws the fly. By integrating Eq. ([3]) over time, we find 
the accumulated yaw, 

Afa = Q J Aipdt w Oil J Aifj dt. (4) 

We confirm that Eq. ([3]) accounts for the observed yaw 
dynamics of flies by plotting Afa versus Co j Ai/jdt for 
10 turning maneuvers in Fig. [4]G. We find that the yaw 
angles of the turns are in excellent agreement with the 
prediction of Eq. (jlj). Thus, by changing the strength 
and duration of the asymmetry in wing rest angles, flies 
can accurately control their turn angle. 

In all forms of locomotion - aquatic, terrestrial, and 
aerial - animals take advantage of mechanical properties 
of their bodies to simplify the complex actuation nec- 
essary to move [lo|, El, [22|. Here, we find a concrete 
example of how this simplification occurs as fruit flies 
modulate their wing motion to generate maneuvers. The 



mechanical properties of the wing hinge appear to be 
fine-tuned to enable insects to modulate their wing pitch 
through only slight active actuation. The spring-like be- 
havior of the wing hinge also connects the timescale of a 
turning maneuver with the timescale of the wing actua- 
tion. Our model predicts, that flight muscles of flies can 
act over several wing beats to bias the pitch of the wings 
and yet generate the sub wing-beat changes in wing mo- 
tion that aerodynamically induce the maneuver. Finally, 
because animals over a wide range of length scales expe- 
rience similar rotational dynamics [19[ , the simple mech- 
anism used by fruit flies may be quite general and should 
likewise simplify the control of flapping flying machines. 
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